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Abstract 

A generic approach for compiling any classical block compression al- 
gorithm into a quantum block compression algorithm is presented. Using 
this technique, compression asymptoticaly approaching the von Neumann 
entropy of a qubit source can be achieved. The automatically compiled 
algorithms are competitive (in time and space complexity) with hand con- 
structed quantum block compression algorithms. 

1 Context 

We are concerned with quantum compression, which is best understood in the 
context of quantum information theory. Informally, the Schumacher coding 
theorem states that a qubit source can be compressed to its von Neumann 
entropy with high fidelity. This statement is strongly analogous to the Shannon 
coding theorem for classical information theory which states that a bit source 
can be compressed to its Shannon entropy with high probability. Quantum 
compression is therefore concerned with algorithms for achieving compression 
near to the von Neumann entropy of a qubit source in a reasonable amount of 
time and space. 

The main result of this paper is a compilation algorithm which takes as 
input any family of classical block compression (and decompression) algorithm 
as well as a description of the quantum source. The output is a quantum block 
compression and decompresion algorithm with the following properties: 

1. If the classical compression algorithm asymptotically approaches Shan- 
non entropy, then the quantum compression algorithm asymptotically ap- 
proaches the von Neumann entropy. 

2. If the classical compression and decompresion algorithm each take space S 
and time T, then the quantum compression algorithm uses space 0(S logT) 
and time 0(±ST 1+£ ) (for any e). 



"This work was also partially completed at ATnT Shannon Labs. Address: Computer 
Science Department, CMU, 5000 Forbes Avenue, Pittsburgh, PA 15213 



Several hand constructed algorithms for block quantum data compression have 
been proposed. These include quantum arithmetic compression 0] which runs 
in 0(n 3 ) time and a combinatoric index compression algorithm pfl running in 
0(n 3 ) time and 0(y/n) space. In further work, the combinatoric index com- 
pression was parallelized to a circuit of only polylogarithmic depth 

The remainder of this paper will discuss some information theory prelim- 
inaries, present the compiler, compare the performance of this compiler with 
hand crafted approaches, and then close with discussion. 

2 (Quantum) Information Theory 

Suppose that we have a qubit source producing qubits in the state |a») with 
probability pj. Then the density matrix of this source will be given by the 
average outer product: 

P = ^Pi\a%){ai\ ='J2<li\ ( l > i)(^i\ 

i i 

where \<pi) are the eigenvectors of the density matrix and t& is the eigenvalue of 

We wish to compress the tensor product of n qubits each drawn indepen- 
dently from this source so as to achieve a high "fidelity". Fidelity is defined 
as the probability that a measurement can distinguish between the input and 
output states. Let \b r ) be the input state chosen with random bits r, and \b r ) be 
the output of a compression/decompression cycle. Then, the fidelity is defined 
by: 

F = E r (b r \b r ) 

where E r is an expectation with respect to the quantum source as well as with 
respect to the outcome of the quantum measurement. 

The Schumacher coding theorem states that the minimum coding rate of any 
compression scheme with a high expected fidelity is given by the von Neumann 
entropy: 

S{p)=Tr{p\ogp) 

Thus, we are particularly motivated by compression schemes which can (asymp- 
totically) achieve the von Neumann entropy. 

3 Generic quantum block compression 

A generic (classical) block compression algorithm, C : {0, 1}" — > {0, l} m , maps 
n bits to to bits with n larger than to. The corresponding decompression algo- 
rithm, D : {0, l} m — > {0, 1}", will then (hopefully) invert the compression to 
give the original state. We can use the general result that any classical func- 
tion can be implemented |^[ (reasonably efficiently) on a quantum computer. In 
particular, If we let 7i, ...,/„ be the bits of the input string and Cj., ...,C m be 
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the bits of the compressed string then we can compile the compression circuit 
into the reversible circuit which calculates: 

Crev(h, •••> I n , Ai, A m ) = I\, ...,/„, A\ @ C\, A m © C m 

where A\, ...,A m are auxiliary bits. The same trick can be done with the de- 
compression circuit to get: 

Drev{Ai, ...,A n ,C\, C m ) = Ai @ Ii, A n © I n , C\, C m 

Pictorially, our circuits are the following: 

I, 

I 2 



A m - 



These elements will be our two fundamental building blocks in the construc- 
tion of generic quantum block compression. Two challenges exist: 

1- Crev and -Drev both have outputs which are entangled with their input. 
This is unacceptable in any quantum compression scheme. 

2. There are some quantum sources with a low von Neumann entropy but 
a large Shannon entropy. We will need to make a circuit which achieves 
compression near the Shannon entropy to achieve compression near the 
von Neumann entropy. 

These two challenges are met in the next subsections. 
3.1 Entanglement Removal 

We are not finished with the compression because the output of Crev leaves the 
compressed bits entangled with the input bits. There is a simple construction 
which removes this entanglement: compose the decompressor with the compres- 
sor to get the following circuit: 



D 



-C0A, C r 

-q0A m c m - 
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Here we have set the auxiliary bits to 0. 

The output of the first compression block is the compressed bits along with 
the uncompressed bits. The uncompressed bits are then fed into the decompres- 
sion block instead of 0's, resulting in an output consisting of the compressed bits 
plus things of the form Ii@D(C(I) )i. This term is precisely when the compres- 
sion is lossless. By assumption, the compression is lossless with high probability 
implying that we will essentially always measure the vector. At this point, 
we can measure all of the "0" qubits. If any of our measurements are nonzero, 
we have detected a low probability event. The process of this measurement 
will project the wavefunction onto the typical subspace, letting us achieve the 
projection done explicitly in Q for "free" as a byproduct of our compression 
process. 

This construction is also particularly elegant because decompression (with 
entaglement removal) can be accomplished by simply reading the circuit back- 
wards. 

In summary, we now have a technique for constructing a quantum compres- 
sion circuit which outputs the compressed string without extra entanglement. 



3.2 Achieving von Neumann entropy 

The above approach only does classical compression of qubits and so it can only 
achieve the Shannon entropy. In particular, it will fail to compress some highly 
compressible quantum state distributions such as the state which is (|0) + |l)) n 
with probability 1. This state is highly compressible because it has a trivial 
distribution over possibilities and yet no classical algorithm will compress the 
qubits well. In order to accomplish compression of this state, we must add an 
extra rotation to the input and output of the compression and decompression. 
For example, the operation H n (Hadamard applied to each bit), would convert 
the input to |0") which is highly compressible. This extra rotation involves 
0{n) operations since it can be done on each symbol from the quantum source 
individually giving us the following circuit: 
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The interpretation of this circuit is simple: we first rotate every qubit to align 
the eigenbasis and the computational basis. Then, we apply the compression 
circuit using the decompression circuit to remove entanglement with the input 
values. Mathematically, we compute 

D Tev (C Tev (R(h), R(I n ), Oi, ..,0 ro )) 

where R is the rotation into the eigenbasis. Once again, decompression is done 
using the circuit backwards. 

This solution works for any quantum source. Given a description of any 
quantum source, the eigenvalue decomposition theorem says that we can con- 
struct a matrix R which satisfies: 

p = RAR T 

where A is a diagonal matrix of eigenvalues and R is a matrix consisting of 
eigenvectors. The R matrix will rotate the input source so that the computa- 
tional basis and the eigenbasis coincide. This rotation does not affect the von 
Neumann entropy of the source. Now, given that the computational basis and 
eigenbasis coincide, the von Neumann entropy will be exactly the same as the 
Shannon entropy. The compiled circuit will therefore achieve the von Neumann 
limit. To summarize, the compiler is the following: 

Make_quantum_compression_circuit(p, C, D) 

1. Compile the quantum reversible circuits -Drev and CVev from the classical 
circuits D and C. 



2. Compose £>rev and CVev as in section 3.1. 

3. Calculate the eigenvalue matrix, R, given p. 

4. Rotate all input and output qubits of the quantum circuit by R as in 



section 3.2 
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3.3 Computational efficiency 



How efficiently can we compile compression algorithms made for classical (irre- 
versible) computers into reversible quantum mechanical programs? There are 
a very few results here which are well summarized in ||. There are a variety 
of tradeoffs all of which are typically parameterized by the space S and time 
T of the original irreversible algorithm. One of the more concrete results is 
that reversible space is SVev = S'logT and the reversible time is (for any e) 
Trev = iflT^ 6 . 

Applying this result to the common 'bzip2' block compression algorithm 
which runs in space S = n and time T = nlogn will require Srev — nlogn 
space and T veY — ^n(n log n) 1+e time (for any e). If we choose e < 1, This is 
a smaller running time than 0(n 3 ) as used by previous hand crafted quantum 
block compression algorithms. It is easy to imagine that hand compilation of 
irreversible compression/decompression algorithms can yield further functional 
improvements. 



4 Conclusion 

The technique for automatically compiling classical compression algorithms into 
quantum compression algorithms gives us a baseline and general structure for 
comparison with hand crafted quantum compression algorithms. The automatic 
approach is reasonably efficient and can give results better than current hand 
crafted approaches in some cases. Hand crafted algorithms which are function- 
ally more efficient than this baseline are a nontrivial improvement. 



The technique of section |3T| applies to any pair of classical algorithms with 
an inversion property: A2(Ai(x)) = x. This includes encryption and decryption 
or encoding and decoding. For algorithms with this property, we can construct 
quantum analogs with a special property: the "inputs" are disentangled from 
the "outputs". This disentanglement comes at a relatively small price: a factor 
of 2 increase in the circuit size. 

In addition to hand crafting efficient compression algorithms, one other sig- 
nificant problem in quantum data compression exists: removing the requirement 
for full information about the quantum source. We know this is (information 
theoretically) possible j?j although no explicit quantum algorithm yet exists 
that achieves asymptotically efficient compression without knowledge of the 
source. Intuitively, this is possible because the precision of the specification 
of R improves exponentially with increasing (qu)bits while the necessary preci- 
sion required in specifying R scales linearly with n, the number of input qubits. 
Consequently, the portion of the bits which must be allocated to specifying R is 
asymptotically 0. An explicit algorithm for doing universal quantum compres- 
sion would yield a significant improvement in our understanding of quantum 
compression. 
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